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^1 Abstract. We present '^€0(1-0) and '^00(2-1) observations of the SA(rs)bc Seyfert 2 galaxy NGC 3147, obtained with the 

^ . IRAM interferometer at V.'9 x V.'6 and V.'6 x l'.'4 resolutions, respectively. We have also observed the central region of NGC 3 147 

?-j ■ with the IRAM 30m telescope (at resolutions of 22" and 12" for '^CO(l-O) and '^C0(2-l), respectively), in order to obtain 

C/3 complete sampling at low spatial frequencies. These observations have been made in the context of the NUclei of GAIaxies 

(NUGA) project, aimed at the study of the different mechanisms for gas fueling of active galactic nuclei (AGN). A central peak 
seen mainly in '-C0(2-l) and a ring-like structure at r ^ 10" ~2kpc dominate the '"CO maps. In '^CO(I-O) an outer spiral at 
r 20" ~4kpc is also detected, not visible in '-C0(2-l) emission because it falls outside the field-of-view of the primary beam. 
^ ' The average hi/Iio line ratio is ~ 0.7 in temperature units over the region mapped in both lines, consistent with the optically 

I thick emission expected in the nuclei of spiral galaxies. The kinematics of the molecular structures are quite regular, although 

00 ■ there is evidence for local non-circular or streaming motions. We show that the molecular gas distribution is similar but not 

exactly identical to those of star formation tracers, i.e., infrared (Spitzer) and ultraviolet (GALEX) emission. This agreement is 
consistent with a scenario of steady, ongoing formation of stars from the molecular clouds at a rate of ~ 1 Mq yr"' within the 
QQ . innermost 4 kpc in radius. 

' Using a near-infrared (NIR) image obtained with adaptive optics at the Canada-France-Hawaii Telescope (CFHT), we identify a 

OO , weak bar in NGC 3147, which is classified as non-barred galaxy in the optical. We then compute the gravity torques exerted by 

' this stellar bar on the gas. The torque is obtained first at each point in the map, and then azimuthally averaged with a weighting 

determined by the gas surface density traced by the CO emission. We find that the gas inside the inner CO ring is subject to a 
net negative torque and loses angular momentum. This is expected for gas at the ultra-harmonic resonance (UHR), just inside 
' the corotation resonance of the stellar bar. In contrast, the gas outside corotation, in the spiral arms comprising the outer spiral 

JH I structure, suffers positive torques and is driven outwards. We conclude that some molecular gas is presently flowing into the 

central region, since we find negative torques down to the resolution limit of our images. 

Key words, galaxies: individual: NGC 3147 - galaxies: spiral - galaxies: active - galaxies: nuclei - galaxies: ISM - galaxies: 
kinematics and dynamics 



1. Introduction 

The N Uclei of GAIaxies (NUGA) project dGarcfa-Burillo et al. 
2003i) is an IRAM Plateau de Bure Interferometer (PdBI) sur- 
vey of nearby active galaxies to map the distribution and dy- 
namics of molecular gas in the inner 1 kpc at high spatial res- 
olution (~ 0'.'5 - 1", coiTesponding to ~ 50 - 100 pc), and to 
study the mechanisms for gas fueling of low-luminosity active 
galactic nuclei (AGNs). 
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Most galaxies possess central supermassive black holes 
(SMBHs), and gas accretion onto these black holes is the phe- 
nomenon usually invoked to explain nuclear activity in galax- 
ies. However, even if most galaxies host black holes, the ex- 
istence of nuclear activity is far from universal. It is not clear 
whether the main limiting factor is the global gas mass avail- 
able for fueling the AGN or the mechanisms for efficiently re- 
moving the angular momentum of the gas. 

The main and non-trivial problem linked to the fueling of 
AGN is the removal of angular momentum from the disk gas 
(e.g.. Jogee„2006. and references therein), a process which can 



2 



Casasola etal.: NGC 3147 



Table 1. Fundamental parameters for NGC 3147. 



Parameter 


Value 


References" 
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lU ID JJ.O 
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X Hubble Type 


J.OO 




Inclination 


zyo 




Position Angle 


150° 


(2) 


Distance 


40.9 Mpc (1" = 198 pc) 


(1) 


Lb 


4.3 X lO'o Lo 


(3) 


Mm 


8.3 X lO'* Mo 


(3) 


Mdust(60 and lOO^m) 


1.7 X 10^ Mo 


(3) 


^FIR 


3.4 X 10'° Lo 


(4) 



(1) NASA/IPAC Extragalactic Dat abase (NED) ; (2) Lyon 
Extragalactic Database (LEDA); (3) ISettoni et ai] j2003l) : (4) 
IRAS Catalog. 



be accomplished through non-asymmetric perturbations. These 
can be perturbations of external o rigin, such as galaxy colli- 

ietal Jl986h . or of 



sions, mergers, and mass accretion (iHeckman 
internal origin due to density w aves, such as spirals or bars, and 
their gravity torques (e.g. lSakamoto et aljl999 : Combesl2001l) . 
In addition to primary bars, fueling processes can be associ- 
ated with more localized phenomena, such as nested nuclear 
bars (e.g. lFriedli & Mar tinet 199 3|), lopsidedness or m - 1 per- 
turbations (e.g. lShu et al.M l990: Garcia-Burillo et al.l l2000l) . or 
warped nuclear disks (e.g.|Schinnerer et al. 2000a b). 

Since molecular gas is the predominant phase of the inter- 
stellar medium (ISM) in the inner kiloparsec of spiral galax- 
ies, the study of its morphology and dynamics represents an 
optimal possibility for investigating AGN fueling mechanisms 
and their link with circumnuclear star formation. In order to 
do this, high-resolution maps of molecular gas are needed. 
Previous CO single dish dHeckman et al. _ 19891: Young et al. 
1995': 'B raine et alJll993l: ICasoh et al jll996l: IV ila-Vilaro et al. 
1998) and interferometr ic^Sakamoto et al.lll9 99: Regan et al. 
200lHHelfer etalJ2003h surveys have mapped the gas in galax- 
ies with relatively low spatial resolution (4" - 7"). Moreover, 
the majority of these surveys only included small numbers of 
AGN in their samples. This paper, dedicated to the galaxy 
NGC 3147, is the tenth of a series where results obtained for 
the galaxies in the NUGA sample are described on a case-by- 
case basis. 

NGC 3147 {D = 40.9 Mp c, \" ^ 198 pc Hq = 
75 km s ' Mpc~') is an iso lated jBettoni et alj l2003b Seyfert 



2 galaxy dHo et al.l Il997h of Hubble morpho l ogical type 



SA(rs)bc. Both ROSAT ("Roberts & Warwick' '2000) and 
Chandra (iTerashkna & Wilso n 2003) observations have shown 
that this galaxy possesses a pointUke nuclear X-ray source, 
which is presumed to be the location of the AGN. The 0.3- 
lOkeV Chandra image reveals a bright, compact source sur- 
rounded by very faint, soft, and diffuse emission. The 2-10 keV 
core is clearly detected and confined to a region of 2" in ra- 
dius. Interferometric observations of the continuum emission 



at cm wavelengths with MERLIN and the VLBA have also 
shown a pointlike nonthermal continuum source at the posi- 
tion of the nucleus of NGC 3147, coincident with the X-ray 
source dUlvestad & H0II2OOII: iKrips et alJEoOel l2007a l). Very 
recent optical and X-ray observations of the nuclear region of 
this galaxy suggest that NGC 3 147 is the first "true" Seyfert 
2 in the sense that it intrinsically lacks a broad-line region 
( Bianchi et al. 20081) . Table [T] reports the fundamental charac- 
teristics of NGC 3 147. 

The structure of this paper is as follows. In Sect. 2, we 
describe our new observations of NGC 3 147 and the litera- 
ture data with which we compare them. In Sects. 3 and 4, we 
present the observational results, both single dish and interfero- 
metric, describing the principal properties of the molecular gas 
including its morphology, its kinematics, and its excitation. A 
comparison between the CO observations and those obtained 
at other wavelengths is given in Sect. 5. In Sect. 6, we describe 
the computation of the gravity torques exerted on the gas by a 
weak stellar bar that we have identified using a near-infrared 
(NIR) image. In Sect. 7, we discuss our principal observational 
and theoretical results, which are summarized in Sect. 8. 



2. Observations 

2.1. IRAM interferometric CO observations 

We observed the 7 = 1 - and 7 = 2 - 1 lines of '^CO in 
NGC 3 147 using the IRAM PdBI in October 2004, with the 
array deployed in its ABCD configurations. The six 15 m an- 
tennae were equipped with dual-band SIS receivers yielding 
SSB receiver temperatures between 40 and 50 K at both fre- 
quencies. The precipitable water vapor ranged from 4 to 10 mm 
(i.e., giving opacities of ~ 0.2 - 0.3), resulting in system tem- 
peratures of approximately 200 - 300 K on average. The spec- 
tral correlators were centered at 1 14. 197 GHz (2793.75 km s"') 
and 228.390GHz (2793.27 km s"'), respectively. The coordi- 
nates of the PdBI phase tracking center are given in Table [T] 
and correspond to the (6 cm VLBA) nuclear radio position of 
Nagaret aLl(l2002h . 



Data cubes with 256 x 256 pixels (0'.' 46 pixel for 
'^00(1-0) and 0'.'23pixer' for i2cO(2-l)) were created over 
a velocity interval of -240 kms"' to -1-240 kms"' in bins of 
10 kms '. The images wer e reconstructed using the standard 
IRAM/GILDASQ software dGuilloteau & Lucasll2000l) and re- 
stored with Gaussian beams of dimensions 1'.'8 x 1'.'6 (PA 
= 62°) at 115GHz and 1'.'4 x 1'.'2 (PA = 61°) at 230GHz. 
In the cleaned maps, the rms levels are 2mJybeam ' and 
4mJybeam-' for the i2cO(l-0) and '^00(2-1) lines, re- 
spectively (at a velocity resolution of 10 kms"'). The con- 
version factors between intensity and brightness temperature 
are 32K(Jybeam-i)-i at 115 GHz and 14K(Jybeam-i)-i at 
230 GHz. 

A continuum point source is detected at both 3 mm and 
1.3 mm. The continuum flux density is 5.2 + 0.2 mJy at 3 mm 
and 2.2 + 0.4 mJy at 1.3 mm. Both values are consistent with 
the synchrotron source detected at centimeter wavelengths for 
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Fig. 1. '^CO(l-O) velocity channel maps observed with the IRAM PdBI+30m in the nucleus of NGC 3147, with a spatial res- 
olution of 1'.'9 X 1'.'6 (HPBW). The center of the maps, given in Table [1] is ajooo = lOMe^SB.e', ^2000 = 73°24'03". Velocity 
channels range from AV = -230 kms"' to +220 kms"' in steps of lOkms"' relative to Vhei = 2813kms"'. The contours begin 
at 5 mJy beam" ' , their spacing is 5 mJy beam" ' , and their maximum is 50 mJy beam" ' . 



a power-law spectrum with a slope of -1. All maps presented 
in the paper are continuum-subtracted. 

2.2. IRAM single dish CO and HCN observations 

We performed IRAM 30 m telescope observations in a 5 x 
5 raster pattern with 7" spacing in July 2002 and in June 
2004. We used 4 SIS receivers to observe simultaneously at 
the frequencies of the '^COCI-O) (115 GHz), the i2cO(2-l) 
(230 GHz), and the HCN(l-O) (89 GHz) lines. The half power 
beam widths are 22" for i2cO(l-0), 12" for '^00(2-1), and 
29" for HCN(l-O). The typical system temperatures were ~ 
250 K at 115 GHz, ~ 350 - 700 K at 230 GHz, and ~ 120 K 



at 89 GHz. The line intensity scale throughout this paper is ex- 
pressed in units of T^^, the beam-averaged radiation tempera- 
ture. The value of Tmb is related to T^, the equivalent antenna 
temperature (corrected for rear spillover and ohmic losses) re- 
ported above the atmosphere, by 77 = T^T^h where 77 is the 
telescope main-beam efficiency, rj is 0.79 at 115 GHz, 0.57 at 
230 GHz, and 0.82 at 89 GHz. All observations were performed 
in "wobbler-switching" mode, with a minimum phase time for 
spectral line observations of 2 s and a maximum beam throw of 
240". The pointing accuracy was ~ 3" rms. 

Single dish observations were used to compute short- 
spacings and complete the interferometric measurements (e.g. 
iCombes et al..2004.) . In particular, short-spacing visibilities are 
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Fig. 2. Same as Fig.[T]but for the '^CO(2-l) line, with a spatial resolution of 1'.'6 x 1'.'4. The velocity range is the same as for Fig 
[1] The contours begin at lOmJy beam their spacing is 5 mJy beam and their maximum is 45 mJy beam"^ 



computed from a map built by interpolation of the 30 m beam 
and multiplied by the PdBI primary beam. We combined 30 m 
and PdBI data, using the SHORT-SPACE task in the GILDAS 
software. To find the best compromise between good angular 
resolution and complete restoration of the missing extended 
flux, we varied the weights attached to the 30 m and PdBI 
data. After writing the combined datasets to visibility tables, 
converting to maps using standard data reduction procedures, 
and deconvolving using the Clark algorithm, we obtained maps 
with angular resolutions of I'.'SS x r.'63 at PA 61.5° for the 
'2CO(1-0) map and 1'.'61 x 1'.'4 at PA 65.4° for the '2cO(2-l) 
map. The weights were adjusted in order to obtain the same 
mean weights in the single-dish data as in the interferometric 
data in the uv range of 1.25 Dl A to 2.5 Dj A (D = 15 m). All 



figures presented in this paper are made with short-spacing- 
corrected data. 

Figs. [Hand |2] display the channel maps for the '^CO(l-O) 
and '^CO(2-l) lines, respectively. All maps are centered on the 
position of Table [1] and the dynamical center coincides with 
this center, which is also the position of the AGN (the radio 
continuum source). 

Figs.|4]and|5]display the single dish data, for the '^CO(l-O), 
'2CO(2-l), and HCN(l-O) lines. The two CO lines have been 
mapped on a 5 x 5 grid with 7" spacings, while the HCN(1- 
0) line has been mapped on a 3 x 3 grid with 7" spacings and 
the nine HCN( 1 -0) spectra have been averaged to improve the 
signal-to-noise. 

To estimate the flux filtered out by the interferomet- 
ric observations, we have computed the flux measured by 
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Fig. 3. Radial brightness profiles are shown in the top panel, 
with J given as a dotted line (blue), and IRAC 3.6 jim as a solid 
one (red). The vertical dashed line corresponds to 200 pc. The 
lower panels display the radial runs of ellipticity, e, position 
angle 0, and cos(40) terms in the ellipse fitting residuals. 



Young et al.l ( Il995h with the FCRAO 14 m telescope (HPBW 
= 45"). They found a '^CO(l-O) intensity towards the cen- 
ter of /(CO) = 5.75 Kkms ' (in the T* scale), which cor- 
responds to an integrated flux of 5 (CO) - 242Jykms ' if 
we adopt a conversion factor of 42JyK"'. In the same re- 
gion, we measure 130Jykms ' in '^CO(l-O) with the PdBI 
alone, and 290Jykms ' with the IRAM 30 m in a beam of 
22". With short-spacings added, we estimate the total flux to 
be 295Jykms ' in '^CO(l-O), which is compatible with the 
FCRAO value within the uncertainties. In '^00(2-1), the PdBI 
alone indicates 30Jykms"' with a large region with negative 
flux levels. With the IRAM 30m in '^CO(2-l) we measure a 
total flux of 140Jykms ' in a beam of 12", and when short- 
spacings are included we recover the value of 1 5 1 Jy km s" ' . 

The implication is that much of the '^CO(2-l) emission 
comes from smoothly distributed gas. 

2.3. Near-infrared observations 

For modeling the gravitational potential, we used NIR images 
which have the advantage of being relatively free of dust ex- 
tinction and more representative of old stellar populations than 
blue or visible images. 

The NIR images were obtained in December 1998 at the 
3.6 meter Canada-France-Hawaii Telescope (CFHT), using the 
CFHT Adaptive Optics Bonnette (AOB) and the KIR infrared 
camera. The AOB, also called PUEO after the sharp-visioned 
Hawaiian owl, delivers essentially difiraction-limited images 
at H and K' (O'.'ll and 0"14, respectively), and images with 
FWHM ~ O'.'l at 7 with guide stars as faint as R = 14 (see 



Rigaut et al.l Il998l) . Here the Seyfert nucleus of the galaxy 
was used for wave front sensing. The KIR 1024 x 1024 pixel 
HgCdTe array has O'.'035/pixel, providing a field of view of 
36" X 36". The observations were obtained in excellent seeing 
conditions (~ 0"6 in the V band), with several images taken in 
a dithering procedure to correct for camera defects. The total 
on-source integration time was 8 minutes each for J and K', 
and the total observing time for the two filters was 50 minutes 
including sky measurements and overheads. 

A J - K map was constructed by first subtracting sky emis- 
sion, as estimated from the outer regions of the images. Then 
mean instrumental magnitude zero points were used to calcu- 
late a J - K image, after aligning them to a common center. 

2.4. Archivai observations witii Spitzer and GALEX 

We acquired public images of NGC 3147 at other wavelengths 
to compare our molecular gas observations with other star for- 
mation tracers. 

We used a far-ultraviolet (FUV) image from the GALEX 
satellite, whose band is centered at /Igff = 1516 A. This im- 
age has been already used and studie d in the context of the 
GALEX Nearby Gala xies Survey (NGS. ICil de Paz et al.ll2004 : 
Bianchiet al]|2003allbl) . a project that spans a large range of 
physical properties and morphological types of nearby galax- 
ies. The image has been obtained with a total exposure time 
of 1693 s and covers a square region on the sky of size ~ 
5500" X 5500", i.e., much larger than the extent of the opti- 
cal disk of NGC 3 147, with r.'5 pixels. As the image was re- 
duced with the GALEX data pipeline, it is akeady expressed 
in intensity units and sky-subtracted. The total FUV calibrated 
magnitude is 14.99 + 0.01, corresponding to a FUV flux density 
of3669 + 17yuJy. 

We also acquired infrared (IR) images obtained with the 
IRAC camera on Spitzer, available thanks to the project 'A 
Mid-IR Hubbl e Atlas of Galaxi es' (Principal Investigator: G. 
Fazio, see also Pahre et al .|2004 . These images range in wave- 
length from 3.6 fim to 8 fim, and were reduced with the Spitzer 
data pipeline (version S 14.0.0). As in the case of the GALEX 
image, they cover a large sky area (~ 1400" x 700"). These 
images have been rotated to canonical orientation (north up, 
east left) in order to allow correct superposition on the molec- 
ular gas map. The pixel size is 1"20 in the final images. The 
stellar component at 3.6yum and 4.5 //m was subt racted from 



the 8.0 jim image ac cord i ng to t he prescriptions of iHelou et al 



(|2004 and Pahre et al.l (|2004|) . This procedure provides a 
"dust-only" image at 8 fim, which we will discuss in Sect. 15.21 

Using the 3.6 //m image, we calculated radial surface- 
brightness profiles by fitting ellipses (using the IRAF/STSDAS 
task ellipse). The center was fixed, but the ellipticity and 
the position angle were allowed to vary. This brightness 
profile is shown in Fig. |3] together with the /-band ra- 
dial profile obtained in a similar way from the CFHT im- 
age. The zero point for the IRAC image was taken from 
the IRAC documentation at the Spitzer Science Center 
(http://ssc.spitzer.caltech.edu/irac/calib/) . 
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Fig. 4. Spectra maps of NGC3147 made with the IRAM 30m with 7" spacing in '2cO(l-0) (top) and '2cO(2-l) (bottom). The 
positions are offsets relative to the center of NGC 3147, whose coordinates are listed in Table [T] Each spectrum has a velocity 
scale from -500 to 500 kms ', and an antenna temperature scale (T^) from -0.05 to 0.25K for '^CO(l-O) and to 0.20K for 
i2C0(2-l). 



No corrections were made to the "standard" IRAC photo- 
metric calibrations. 

Also shown are the runs of ellipticity, ellipse position an- 
gle, and cos(40) residuals of the ellipse fitting. A local peak in 
the ellipticity profile at R ^ 5" (1 kpc) at constant position an- 



gle (85°) corresponds to a weak bar/oval feature that will be 
discussed in Sect. 15. II 



Casasolaetal.: NGC 3147 



7 



TAt (mK) 



n 



Velocity (km/s) 

Fig.5. HCN(l-O) spectrum towards the center of NGC3147, 
averaged over the 9-point map made with the IRAM 30 m with 
7" spacing. 



3. Single dish results 

The observations performed with the A and B receivers of the 
IRAM 30 m telescope in the two '^CO lines covered the central 
~50", corresponding to the central ~10 kpc (in diameter) of the 
galaxy (Fig.|4]i. The 25 observed positions show that the central 
region of NGC 3147 presents extended molecular emission in 
both '2cO( 1 -0) and ^^CO{2- 1 ) (Fig.gli. The maximum detected 
T„,b is 0.32 K in '^COCI-O) and 0.31 K in i2cO(2-l), in the 
southeast region corresponding to the offsets (10", -10") and 
(5", -10") relative to the galaxy center 

The combination of observations of the '^CO(l-O) and 
'^CO(2-l) lines allowed us to compute the line ratio, /?2i = 
hi/Iw = / T,„i,i2 - l)dv/j T,„hil - 0)dv, an indicator of phys- 
ical properties of the molecular gas such as excitation and opti- 
cal depth. After smoothing the '^CO(2-l) data to the '^€0(1-0) 
beam resolution (22") and correcting for the different beam ef- 
ficiencies, we found a mean line ratio for the central region of 
NGC 3 147 of 0.8. This value is consistent with the optically 
thick emission from mo lecular gas expected for the central re- 
gions of spiral galaxies ( Braine & Combeslll992 ). 

Assuming a CO-H2 conver sion factor X - N(H^)/Ic n = 
2.2 X 10^" cm -2 (K km s ')-' jSolomon & Barrettlll99lb . the 
total H2 mass derived from these single dish observations is 
- 8.2 X 10'^ M0, within a radius of 5 kpc. 

The HCN(l-O) line has been observed for 9 positions with 
7" spacing covering the central 43", corresponding to the cen- 
tral ~8.5 kpc (in diameter). The HCN(l-O) emission has been 
detected most significantly in the northwest part of the ob- 
served grid corresponding to the offset (-7", 7"). The HCN(1- 
0) average spectrum over the 3x3 grid is plotted in Fig. |5] The 
CO(1-0)/HCN(1-0) ratio is equal to 20 on average. Because the 
observed region contains the AGN, we would expect enhanced 
HCN emission there. Hence, the observed CO(1-0)/HCN(1-0) 
ratio is unexpectedly high. For instance in NGC 1 097, this ratio 



is 3 in the nucleus and 10 in the star-forming ring ( Kohn o et al 



l2003i) . NGC 3 147 is more similar to the starburst in NGC 695 1 



dKrips et alj2007bh . where the ratio is 30 in the ring, and 2.5 in 
the nucleus. 

To study the physics and distribution of the molecular gas in 
detail we need to analyze both single dish and interferometric 
observations. The analysis in the remainder of the paper will be 
performed on the combined IRAM PdBI-H30m datasets. 

4. Interferometric results: Molecular gas 
properties 

4.1. Morphology and mass of the CO rings 

The '^CO(l-O) and ^^CO(2-l) integrated intensity distiibutions 
are shown in Fig. |6] The figure reveals a nuclear peak, a sym- 
metric and complete inner ring at a distance of about 10" from 
the nucleus, and a larger and incomplete ring at a distance of 
about 20" from the nucleus (left panel). In '^CO(2-l) the nu- 
clear peak is stronger and more visible (right panel), but there 
is no outer ring at a radius of 20" because of the restricted field- 
of-view (FWHM of primary beam of 22"). The outer ring-like 
structure can also be interpreted as spiral arms, more tightly 
wound in the inner regions than in the outer ones. The ra- 
dial distribution of the azimuthally averaged '^CO(l-O) and 
'^CO(2-l) intensities is shown in Fig. |7] where there are clear 
dips between the central peak and the inner ring and between 
the inner ring and the outer spiral structure (partial outer ring). 

Both '^CO(l-O) structures are clumpy, composed of indi- 
vidual giant molecular cloud complexes, each with a mass of 
a few 10^ to 10** M©. The total H2 mass derived from the in- 
terferometer maps, assuming the same CO-to-Hi conversion 
value used before, is M//, = 3.8 x 10** Mq. The inner ring alone 
contributes more than half of the total masses; M//, = 2 x lO' 
Mq, and the two parts of the outer spiral have similar mass, 
the western part has a mass of M^^ - 9.5 x 10^ Mq and the 
eastern one Mh, = 8.5 x 10^ Mq. The total M^, (= 3.8 x 10'' 
Mp)) we find is in good agreement with the mass derived by 
Young et al.l (Il995h for a position centered on the galaxy nu- 
cleus, Mh2 = 4.4 X 10^ Mq. Our mass estimate from single dish 
observations presented in Sect.[3]is larger, since it corresponds 
to a more extended region. In any case, this molecular gas mass 
is very large, more than three times larger than the most mas- 
sive of the othe r NUGA g alaxies studied so far (NGC 2782 and 
NGC 4569, see lHunt et al. 2008; Boone et al. 2007). 

The comparison between the '^CO maps of the two tran- 
sitions, at the same resolution and with the same spatial fre- 
quency sampling, gives insight about the excitation conditions 
of the molecular gas locally, pixel by pixel. Fig. [8] shows that 
when the '^CO(2-l) data are tapered and convolved to the 
'^CO(l-O) resolution, the maxima of the '^CO(2-l) ring agree 
quite well with the '^CO(l-O) peaks. The ratio map is shown 
in Fig. |9] where the values range from 0.1 to 1.8. The mean 
ratio is ~0.7, consistent with the value found with single dish 
observations and also, in this case, consistent with the bulk of 
emission being optically thick, as expected in the nuclei of spi- 
ral galaxies. Some regions, especially the center, reach a higher 
ratio (~ 1.3- 1.4) where the gas is warm and dense, while in the 
disk and spiral arms it is colder and more diffuse. The peak in 
the northwest region of the CO(2-1)/CO(1-0) ratio map is co- 
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Fig. 6. Left panel: '^CO( 1-0) integrated intensity contours observed with the IRAM PdBI+30 m toward the center of NGC 3 147. 
The cross marks the coordinates of the center as given in Table [T] with offsets in arcseconds. The rms noise level is cr = 
0.15 Jy beam"' km s"'. The map, derived with 2cr clipping, has not been corrected for primary beam attenuation. Contour levels 
run from 4cr to 17cr with 2.5cr spacing. In this map the full +240 km s"' velocity range is used. The beam of l'.'9 x V.'6 is plotted 
at lower left. Right panel: Same for '^CO(2-l). The rms noise level is cr = 0.2 Jy beam"' km s"'. Contour levels run from 5cr to 
18cr with 2.5cr spacing. The beam of I'.'A x I'.'l is plotted at lower left. 
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Fig. 7. Left panel: Radial distribution (azimuthal average, deprojected to face-on orientation) of the '^CO(l-O) integrated intensity 
shown in Fig.|6] Right panel: Same for '^CO(2-l). 



incident with the strongest HCN(l-O) emission corresponding 
to the offset (-7",7"). 



4.2. Kinematics 

The channel maps in Fig. [T| display overall a gross regularity 
of the large scale kinematics, following the expected spider di- 
agram. However, there are some local wiggles superposed on 
this pattern, i.e., streaming motions. Such streaming motions 
are also visible in Fig.[TOl which shows the isovelocity curves 
(first-moment map) of the '^CO(l-O) emission superposed on 
the '^CO(l-O) integrated intensity. Some of the perturbations 



in the velocity field are clearly coincident with the spiral arms, 
as expected from density wave theory. 

A position-velocity (p-v) cut along the major axis using a 
position angle of 150° is presented in Fig.[TT] The major-axis 
p-v diagram shows approximately regular kinematics, although 
there are clear signs of deviations from "flatness" at ~15"on 
either side of the nucleus of amplitude ~40kms"'. 

The minor-axis p-v diagram is displayed in Fig. [12] The 
wiggles are more conspicuous on the major axis in Fig.fTTIthan 
in Fig. [12] Although the minor-axis velocities are not com- 
pletely symmetric about the nucleus, they hint that the outer 
spiral arm/pseudo-ring at a radius of 20" is already outside 
corotation. This is because the sign of the radial streaming 



CasasolaetaL: NGC 3147 



9 




Fig. 8. i2CO(l-0) contours as in Fig. |6] (left panel) superposed 
on the greyscale '^CO(2-l) map that has been tapered, con- 
volved to the same resolution, and corrected for primary beam 
attenuation, in units of Jybeam"' kms The beam is plotted 
at lower left. 
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Fig. 9. Contours and greyscale of the CO(2-1)/CO(1-0) ratio 
map. Contours run from 0. 1 to 1 .8 in steps of 0. 1 in temperature 
units. 



motions (visible on the minor-axis) is expected to change at 
corotation, while the tangential streamin g (visible on the major- 
axis) does not change at corotation (e.g. Icanzian|[l993l) . Thus, 
we conclude (although with substantial uncertainty) that coro- 
tation occurs at ~15" (3kpc), consistent with the minor-axis 
p-v plot. This result will help constrain the pattern speed of the 
weak bar/oval which we discuss in Sect. [T] 
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Fig. 10. Overlay of the integrated '^CO(l-O) emission in grey 
scale, same as Fig. |6] (left panel), with the CO mean-velocity 
field in contours spanning the range -200 km s ' to 200 km s ' 
in steps of 10 km s"' . The velocities are referred to V/,(.; = 2813 
km s '. Solid (dashed) lines are used for positive (negative) 
velocities. The dotted line indicates the major axis (PA = 150°). 

4.3. Dynamical mass 

The peak velocity of -200 kms ' in the major-axis p-v dia- 
gram occurs at a radius of ~7-8" (1.5 kpc). This value, when 
corrected for disk inclination (we adopt / = 29.5°, see Table 
[U by a factor sin /, becomes 400 km s"'. From the corrected 
peak velocity, we can estimate the dynamical mass within this 
radius: M(R) = a x where M(R) is in Mq, R in kpc, and 

V in km s"' . Assuming the most flattened disk-like distribution 
(a — 0.6), gives a dynamical mass Mdyn = 3.3x10'° Mo within 
a diameter of ~3 kpc. Assuming a roughly flat rotation curve 
(see Sect. this would imply a mass Mdyn = 8.9 X 10'° Mo 
within 4 kpc (roughly the radius of the PdBI primary beam for 
our '^CO(l-O) observations). With a total molecular gas mass 
of - 4xlO''Mo(seeSect.|4l]i, we would thus estimate a molec- 
ular mass fraction of 4 - 5%. 



5. Comparison with other wavelengths 

It is instructive to compare the distribution of the molecular 
gas in the inner kpc of NGC 3147 with observations at other 
wavelengths. In particular, molecular clouds are thought to be 
the birth site of future generations of stars, and thus should 
be strongly connected to other emission tracing star formation, 
such as far-ultraviolet and warm dust emission. 



5.1. The bar and the lopsided nucleus 

The CO contours are superposed on the CFHT 7-band image 
in Fig. [13] There is a weak bar/oval feature at a position angle 
PA - 85°, just contained inside the inner CO ring. Only the 
NIR reveals this bar, since NGC 3147 is not classified optically 
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as a barred galaxy. Its late type (Sbc) corresponds to a non- 
dominating bulge, and the bar feature is likely not diluted by 
this bulge. The major axis orientation differs by 65° from the 
bar PA, which minimizes confusion due to projection effects. 



Mollenhoff &Heidtl ( 1200 Ih have obtained JHK images of 



NGC 3147 with the 2.2 m telescope of the Calar Alto observa- 
tory using the MAGIC NICMOS3 camera, together with 40 
other unbarred galaxies with low inclinations. They decom- 
posed the images into bulge and disk components, and noticed 
that the bulge in NGC 3147 appears at a very different posi- 
tion angle than the disk; they concluded that the bulge must be 
triaxial, which was a rare feature in their sample. We have re- 
reduced their MAGIC J and K' images and the oval is present, 
similar to the feature seen in our CFHT images, although with 
noisier isophotes. Our kinematic data suggest that this bar may 
be the agent of the inflow that we discuss in Sect.|6] 
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Fig. 11. '^CO(l-O) position-velocity (p-v) diagram along the 
major axis of NGC 3147 using the whole velocity range +1- 
240 km s ^ The coordinates of the center (0,0), the heliocentric 
velocity, and PA used are given in Table [T] Contour levels are 
from 3.2 to 32mJy beam ' in steps of 3.2mJy beam 



The ellipticity profile shown in Fig. [3] (see Sect. I2.4l i helps 
us estimate the length of the bar/oval distortion. The ellipticity 
peaks at ~5" (1 kpc), but the bar extends roughly to a radius 
of ~7.5" (1.5 kpc). This is also seen in the 7-band image (Fig. 
\T3[ where the bar is roughly 16" in diameter, and in the Fourier 
decomposition of the stellar potential discussed in Sect. |6] We 
therefore adopt a bar radius of ~7.5" or ~ 1 .5 kpc, although this 
estimate is subject to a rather large uncertainty, +30%. 

The J - K image of NGC 3147 is shown in Fig. [H with 
the '^CO(2-l) integrated intensity overlaid as contours (see 
Sect. 14. 11 1. Coincident with the nuclear CO emission, there is a 
clearly-defined non-axisymmetric structure of ~700 pc in size, 
which is bluer in J - K than the surrounding regions that are 
devoid of gas. This central asymmetric "blob" or "lopsided nu- 
cleus" could be the destination of the gas inflow that we discuss 
in Sect.m 



5.2. The inner ring and tlie outer spiral 

Fig. [15] displays the FUV image obtained with the GALEX 
satellite, described in Sect. 12.41 Inspection of Fig.fTSlshows that 
NGC 3147 has a quite well-defined outer spiral structure, out- 
side the range of our CO observations. There is also a compact 
nucleus with a central peak, surrounded by a roughly axisym- 
metric "plateau" of emission which breaks up into a pseudo 
ring at a radius of ~20". 

Fig. [16] shows the superposition of GALEX FUV image 
with the '2CO(1-0) (left panel) and '^00(2-1) (right) integrated 
intensity contours. The cross shown in Fig.[T5]is at the center of 
our CO observations, and coincides with the peak of the FUV 
nuclear emission. This agreement suggests that the FUV and 
CO emission are cospatial in the nuclear region. 
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Fig. 12. '^CO(l-O) position-velocity (p-v) diagram along the 
minor axis of NGC 3147 using the same velocity range, center, 
PA, and systemic velocity as for Fig. [TT] Contour levels are 
from 2.8 to 28 mJybeam"' in steps of 2.8 mJy beam"'. 

However, there is no clear ring in the FUV image. The inner 
1 0"-radius CO ring is contained inside the circular region of the 
FUV emission, and the outer CO spiral (20"-radius), at least its 
eastern part, has some overlap with the FUV spiral, but there 
are clear shifts between the two distributions. The '^CO(2-l) 
emission (essentially the inner ring) also appears to settle inside 
the GALEX FUV inner "plateau" (Fig. [IS]). 

The "outer spiral" in '^CO(l-O) falls in the interarm 
"gap" in the GALEX image. Offsets between FUV and 
CO (or FIR and HI) emissions are not u ncommon in 
grand design spiral galaxies, su ch as M 100 dRandl 1995 : 
Sem pere & Garcfa-Buriliol Il997h and M 51 dCalzetti et al 
.2005) , and may relate to star-formation efficiency and 
timescale variations in response to a spiral density wave. The 
star formation and all related tracers are often located in differ- 
ent regions of the spiral arms: FUV emission is more promi- 
nent at the outer edge of the arms, where the dust extinction is 
low, while FIR and Ha emission is stronger at the inner edge. 
However, in NGC 3147 we find CO in the middle of the inter- 
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Fig. 13. Left panel: '^CO(l-O) contours (0.6 to 2.6 by 0.35 Jy km s ' beam ') superposed on the near-infrared J image from the 
CFHT, shown in logarithmic levels. The central 50" are shown. Right panel: Same for '^CO(2-l) contours (1.0 to 3.5 by 0.5 Jy 
km s ' beam '). The central 30" are shown. 



arm region in the GALEX image, not at the inner or outer edge 
of the arm; this result is quite rare. 

Unlike the FUV, the 8 fim "dust-only" image described in 
Sect. l2.4l shows an inner ring apparently associated with molec- 
ular gas. 

Fig. [17] shows this image, which clearly exhibits an inner- 
ring structure at about 10" radius, together with a central peak. 
Such rings are frequently observed in barred galaxies with 



Spitzer (e.g. iRegan et al.l 12006) . Fig. [18] shows the dust-only 



S/im image with the '2cO(l-0) (left panel) and '^00(2-1) 
(right) intensity contours superposed. The ring visible at 8 fim 
corresponds to both the '^CO(l-O) inner ring and the '^€0(2-1) 
one. However, the outer '^CO(l-O) spiral is slightly inside the 
corresponding structure in dust emission. This may be a simi- 
lar phenomenon to that seen in the GALEX FUV image, where 
molecular gas peaks in the interarm regions as described above. 

We can use the 8 pm dust-only image, together with the 
CO emission itself, to define the size and width of the inner 
gas (and dust) ring. The ring is rather broad, extending over 
roughly 2-3" in radius, or ~500pc. Measured to the middle of 
its width, the inner gas ring has a radius of ~ 9 - 10" (1.8- 
2 kpc) , although it is slightly elongated, being longer along the 
horizontal axis. If instead we measure to the innermost portion 
of the ring, we find a radius of ~ 7" - 8" (~ 1.5 kpc). This 
uncertainty (20-30%) will influence our discussion in Sect. [7] 
where we will use the size of the ring in the context of the 
rotation curve to help constrain the pattern speed and location 
of the resonances that define the kinematics in NGC 3147. 

5.3. Star formation 

The agreement, where present, between the molecular gas dis- 
tribution and the ultraviolet and infrared emission is consistent 
with on-going star formation within the molecular clouds. The 
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Fig. 14. TheJ-K image from the CFHT with '^00(2-1) in- 
tegrated intensity overlaid in contours (1.0 to 3.5 by 0.5 Jy km 
s"' beam"'). North is up and east to the left. 

star formation rate (SFR) can be computed f rom both UV and 
FIR emission using the calibrations given by [KennicuttI ([ 1 998h 
for a 0.1 - 100 Mq Salpeter initial mass function: 



SFR[Moyr-i] = 1.4 x 10"^** Lj/v 

where Luv is the UV luminosity in ergs s" 
wavelength range 1500-2800 A, and 



SFR [Moyr-i] = 4.5 x 10"'*^Lf« 



(1) 



Hz over the 



(2) 
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Fig. 15. GALEX FUV sky subtracted image of NGC 3147. The 
central 120" of the galaxy are shown. The box shows the cen- 
tral 50", field of view of our '^CO observations. 
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Fig. 17. Spitzer/IRAC "dust-only" image at 8 yum of NGC 3147. 
The central 120" of the galaxy are shown. The box shows the 
central 50", field of view of our '^CO observations. 

where L/t/r, expressed in ergs s"\ refers to the total infrared 
luminosity integrated over the entire mid and far-lR spectrum 
(8-1000 yum). We first compute the UV luminosity from the 
GALEX map inside the FOV of 42", to better separate the SFR 
in the PdBl FOV from the total, and find a SFR ~1 Mq yr '. 

We then use the Spitzer MIPS images available from the 
archive to estimate how much FIR flux falls within the FOV of 
the PdB primary beam compared to the total integrated over the 
much larger IRAS beams. Rough measurements give a fraction 
of ~12% for 160/vm and ~22% for 24pm (the 70pm image is 
not available). We then estimate the FIR luminosity by scaling 
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Fig. 19. Strengths (Qi, Q2, and total Qj) and phases ((pi and 
<p2) of the m - 1 and m - 2 Fourier components of the stellar 
potential inside a field of 48" in diameter (r < 4.6 kpc). 

the IRAS 100 /zm and 60jum fluxes by 12% and 119^ respec- 
tively. This gives 1.39 Jy at 60 pm and 3.55 Jy at 100 ;um, which 
results in an estimate for the SFR of ~0.8MQyr"', similar to 
that obtained from GALEX. 

The compatibility between the two estimations of SFR 
(from UV and FIR) means that there is not significant obscu- 
ration, which is probably due to the nearly face-on view of the 
galaxy. 

The SFR per unit area inferred from these values is 
O.O16M0yr ' kpc"^, high for spiral disks but similar to, al- 
though slightly lower t han, typical valu es for IR-selected cir- 
cumnuclear starbursts ( Kennicut3ll998l) . Moreover, given the 
g as mass we derived in Sect. 14.11 and the correlations shown 
in lKennicuttI (1 1998b . the global star-formation efliciency within 
the inner 8 kpc is low, ~l% rather than the more typical value 
of ~10%. On the other hand, the molecular gas is concentrated 
into well-defined structures; more than half the molecular mass 
is found in the inner ring alone (see Sect. 14. 11 1. If we consider 
only the molecular gas in the inner ring, the gas surface den- 
sity is ~320 Mq pc"^, roughly 4 times higher than we obtain by 
averaging over the entire FOV. Nev ertheless, it is sign ificantly 
lower than that in M 5 1 spiral arms dHenrv et alj|2003b . 



- 17% is a rough interpolation of the central concentration we 
would expect at 60yum, given the more compact emission at 24 pm 
and the more extended emission at 160//m. 
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Fig. 16. Left panel: ^^CO(l-O) contours (0.6 to 2.6 by 0.35 Jy km s ' beam superposed on the GALEX FUV image of 
NGC 3147, shown in logarithmic levels. The central 50" are shown. Right panel: Same for i2cO(2-l) contours (1.0 to 3.5 
by 0.5 Jy km s"' beam"'). The central 30" are shown. 




RA offset (arcsec) 



Fig. 18. Left panel: '^CO(l-O) contours (0.6 to 2.6 by 0.35 Jy km s ' beam ') superposed on the 8 ;um "dust-only" (Spitzer/IRAC) 
image of NGC 3147, shown in logarithmic levels. The central 50" are shown. Right panel: Same for '2cO(2-l) contours (1.0 to 
3.5 by 0.5 Jy km s"' beam"'). The central 30" are shown. 



6. Computation of the torques 

6.1. Evaluation of the gravitational potential 

As d escribed in previous papers (e.g., iGarcia-Burillo et aL 
2005h . we assume that NIR images can give us the best approxi- 
mation for the total stellar mass distribution, being less affected 
than optical images by dust extinction or by stellar population 
biases. We here recall the essential definitions and assumptions. 

The y-band image was first deprojected according to the 
angles PA = 150° and ; - 29.5°. We have not deprojected the 
bulge separately, since we do not know its actual flattening, but 



because the galaxy is of late type (Sbc), it does not have a large 
contribution. The image is completed in the vertical dimension 
by assuming an isothermal plane model with a constant scale 
height, equal to ~ 1/1 2th of the radial scale-length of the image. 
The potential is then derived by a Fourier transform method, as- 
suming a constant mass-to-light (M/L) ratio. The M/L value is 
selected to reproduce the observed CO rotation curve. The po- 
tential was also calculated from the lower-resolution MAGIC 
images, which have a much larger field-of-view (162"). The 
same results are obtained, and the potential is shown to be ax- 
isymmetric at large radii. Beyond a radius of 20", the mass 
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density is set to 0, thus suppressing any spurious m - 4 terms; 
this is sufficient to compute the potential over the PdB '^CO(l- 
0) primary beam. 

For the non-axisymmetric part, the potential <!>(/?, 6) is then 
decomposed in the different m-modes; 

CO 

<D(/?, 6) = ^o{R) + Yj W cosCme - 4>„,{R)) 

m=l 

where (b,„iR) and (pmiR) represent the amplitude and phase of 
the w-mode. 

The strength of the m-Fourier component, Qm{R) is defined 
as Qm{R) - ni<^mlR\F Q{R)\, i.e. by the ratio betw een tangential 
and radial forces ("e.g. lCombes & Sanderslll98lh . The strength 
of the total non-axisymmetric perturbation is defined by: 



Qt{R) 



Tpmax 



(R) 



FoiR) 



where Fj"^{R) represents the maximum amplitude of the tan- 
gential force and Fq{R) is the mean axisymmetric radial force. 
Fig. [19] displays these Q values as a function of radius for 
NGC 3147. A bar (m - 2) can be seen clearly, with a peak 
radius of ~ 1 kpc (5"), and a total extent of 1 .5 kpc (7.5") in ra- 
dius. These dimensions are similar to those inferred from the 
ellipticity distribution discussed in Sect. 15.11 and strengthen 
our conclusion that this structure is "bar-like". There is also an 
asymmetry in the m - 2 profile towards the center. The lopsid- 
edness has a correspondence in the gas morphology described 
in Sect. 14. H and the J - K image discussed in Sect. 12.41 The bar 
is regular in phase, but its strength is relatively modest. 

6.2. Evaluation of gravity torques 

After having calculated the forces per unit mass (Fj^ and Fy) 
from the derivatives of <i>(R, ff) on each pixel, the torques per 
unit mass t(x,y) can be computed by: 

t(x,y) = xFy-yF^, 

The torque map is oriented according to the sense of rotation in 
the plane of the galaxy. The combination of the torque map and 
the gas density E map then allows us to derive the net effect on 
the gas at each radius. The gravitational torque map weighted 
by the gas surface density t(x,y) x l,{x,y), normalized to its 
maximum value, is shown in Fi gs . l20l and 1211 for CO(l-O) and 
CO(2-l), respectively. 

To estimate the radial gas flow induced by the torques, we 
first compute the torque per unit mass averaged over azimuth, 
using I,(x,y) as the actual weighting function, i.e.: 



t(R) 



J^nx,y)x(xFy-yF,) 



By definition, t(R) represents the time derivative of the spe- 
cific angular momentum L of the gas averaged azimuthally, i.e., 
t{R)=dL/dt \g. To derive dimensionless quantities, we normal- 
ize this variation of angular momentum per unit time to the 
angular momentum at this radius and to the rotation period. We 
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Fig. 20. CO(l-O) contours overlaid on the gravitational torque 
map (t(x,y) x E(x,y), as defined in text) in the center of 
NGC 3147. The torque map (grey/color scale) is normalized 
to the maximum absolute value of the torques. The derived 
torques change sign as expected in a butterfly diagram, delin- 
eating four quadrants. The orientation of quadrants follow the 
bar orientation in NGC 3147. In this deprojected picture, the 
major axis of the galaxy is oriented parallel to the abscissa. 
The line reproduces the mean orientation of the bar (PA = 20° 
on the deprojected image and 85° on the projected one). 
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Fig. 21. Same as Fig.|20]for the '^CO(2-l) emission taken as a 
tracer of gas surface density. 



then estimate the efficiency of the gas flow as the average frac- 
tion of the gas specific angular momentum transferred in one 
rotation (Trot) by the stellar potential, as a function of radius, 
i.e., by the function AL/L defined as: 



AL 
T 



dL 
dt 



1 

X — 

L 



_ t(R) 

X Ti gt — — X Trot 
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where Lg is assumed to be well represented by its axisymmetric 
estimate, i.e., Lg - R x Vrot- The AL/L curves for NGC 3147 
derived from the CO(l-O) and the CO(2-l) data are displayed 
in Figs.l22l 

The gravitational torque maps in Fig. |20] show that the de- 
rived torques change sign following a characteristic 2D butter- 
fly pattern. The CO contours reveal that for the material in the 
inner ring (radius 2 kpc, 10"), part of the gas is trailing the bar, 
while part of the gas is leading it. The observed gas distribu- 
tion is representative of the time spent by a molecular cloud 
on a typical orbit at this location. The azimuthal average of the 
torques shown in Fig.|22]suggests that the torques are predom- 
inantly negative inside a radius of ~2 kpc, while they become 
zero or positive outside. 

If the near side is the northeast side, the rotation sense of 
the galaxy is counterclockwise, and the whole spiral structure is 
trailing with respect to the rotation. In the spiral structure wind- 
ing into the outer '^CO(l-O) ring, the torques almost cancel out 
on average but are still predominantly positive. Inside the inner 
ring (r < 2 kpc), however, the dominating torques are negative, 
although weak in absolute value. The '^CO(2-l) emission is 
more clumpy, and our map suffers from lower sensitivity; there 
is relatively more emission in the central part. Although the 
emission is spread over three quadrants in the center, the nega- 
tive torques are still predominant there (Fig.l2TI). The negative 
torques there are due to the central asymmetric CO distribu- 
tion, mainly the nuclear component, at radii inside 3-4". The 
component is resolved in '^CO(2-l), and is not aligned with 
the nuclear bar, but rather shifted to the leading side; therefore, 
negative torques dominate. 

7. Discussion 

We find that the two molecular structures (inner ring and outer 
spiral), in addition to a central peak of '^CO(2-l) emission, 
dominate the CO map of NGC 3147's inner region. The inner 
ring coincides with the ring in the dust-only 8 fim image, and 
the lopsided central peak falls at the same position as the asym- 
metric nuclear structure in the J-K map. These features are not 
present in the stars (e.g., J band, IRAC 3.6;um, and GALEX), 
but are rather features of the ISM. On the other hand, the outer 
spiral also appears in the stars as shown by the GALEX image; 
however, there is an offset between them, as the CO emission 
falls in the interarm region relative to the outer spiral (see Fig. 

A probable explanation is that the CO outer spiral is not a 
permanent structure but is transient. The gas participates in a 
spiral wave and forms stars, and then the newly formed stars 
decouple from the gas component. The young stars heat only 
the dust present at their location, while the bulk of the gas has 
moved to another density condensation. The inward migration 
of the dense gas with re spect to the initial resonance has for 
been discussed by, e.g., iRegan & Teuben I (l2003h . This phe- 
nomenon, because it is very transient, is statistically rarely en- 
countered. In this scenario, the CO clouds would not trace stars 
of the same age as those that have heated the dust observable 
at 8 i^m or of those whose emission is detected in the FUV, but 
rather they will be the birth site of future stellar generations. 




R (kpc) 

Fig. 23. Model rotation curve and derived frequencies Q, Q - 
k/2, and Q + k/2, for NGC 3147. The star symbols are the data 
derived from the CO kinematics, deprojected with an inclina- 
tion of /= 29.5°. The horizontal line corresponds to O,, = 125 

km s"' kpc"^ 



For the torques computation, we interpret the results by 
identifying the resonances with the bar Let us note that CO 
resonant rings have already been observed in weakly barred 
galaxies, where the bar is inconspicuous in the optical (e.g. 
NGC 5005, NGC 7217: ISakamoto et al]|2000t ICombes et al 
,2004) . The rotation curve derived from the CO kinematics is 
rather peaked in the center, and corresponds well to the gravi- 
tational model obtained from the NIR image. Both are fitted by 
an axisymmetric mass model, with a rotation curve as shown in 
Fig. |23] which allows us to derive the characteristic dynamical 
frequencies. 

If the stellar bar ends at or just inside its ultra-harmonic 
resonance (UHR), as is canonical for bar dynamics, this means 
that the pattern speed of the bar would be roughly Q.p - 
120 - 130 km s"' kpc"', corresponding to a bar with one or 
two inner Lindblad resonances, as is common. The presence 
of these ILRs is not certain, however, given the possible non- 
circular motions in the bar The inner CO ring would then 
correspond to the UHR (r = 2kpc, 10") just inside corotation, 
which we estimated to be at a radius of ~15" (3 kpc) in Sect. 
14.21 The measurements of the inner ring radius and bar length 
suggest that the bar/oval ends roughly at the onset of the inner 
ring. However, the considerable uncertainties of these measure- 
ments and the finite width of the ring lead to an uncertainty in 
the pattern speed Q.p of the bar 

The outer CO spiral corresponds to the spiral structure 
emerging from the bar, beyond its corotation. The pitch angle 
of these spiral arms is quite small, as is expected for a weak 
bar In this scenario, the gas in the inner CO ring is inside coro- 
tation; due to the negative torques there, the gas loses angular 
momentum and is therefore flowing inwards. Due to the weak 
bar, the rate of inflow is relatively low, with a time-scale of 
about 5 rotations (Fig.l22li. 

The amount of gas flowing inwards is consistent with what 
might be needed to fuel the Seyfert 2 nucleus. The gas surface 
density in the central kiloparsec is about 10 Mopc"^; if the 
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Fig. 22. The torque, or more precisely the fraction of the angular momentum transferred from/to the gas in one rotation -dLfL- 
is plotted for '^COCI-O) (left) and '2CO(2-l) (right). 



gas at 1 kpc from the center falls in within 5 dynamical time- 
scales (of ~20 Myr), then the feeding rate is about 0.5 M© yr 
Assuming a mass-to-energy conversion efficiency of e ~ 10% 
(L = dM/dtc^e), then a luminosity of order 3 xlO"^^ erg s"' is 
accounted for, more than enough for this low-luminosity AGN. 

The spatial resolution of the CO maps is not yet sufficient 
to determine more precisely whether the gas is actually feeding 
the AGN, but we can at least conclude that we are observing the 
gas inflow into the ~1 00-200 pc scale circumnuclear region. 

8. Conclusions 

The molecular gas in the Seyfert 2 galaxy NGC 3147 has been 
mapped with high resolution (1'.'9 x 1'.'6 Gaussian beam for the 
'2CO(1-0) line) inside a radius of 25" (~5 kpc). The CO emis- 
sion shows a central peak (mainly in '^CO(2-l)) and a ring at 
distance of 2 kpc. In '^CO(l-O) also an outer spiral at distance 
of 4 kpc is detected. The observed CO has a mean line intensity 
ratio /21//10 ~0.7, consistent with the optically thick emission 
expected in the nuclei of spiral galaxies, and broadly regular 
kinematics with some evidence for local non-circular motions. 

Comparing the molecular gas distribution with tracers of 
star formation, we find that central emission and the inner CO 
ring coincide well with the FUV (GALEX) and 8 yum (Spitzer) 
emission, while the outer CO spiral structure/ring is instead 
located within an interarm interval detected by GALEX and 
Spitzer. We interpret this partial disagreement as due to the 
transient nature of the outer CO ring. It would be a structure 
in re-condensation that traces a future stellar generation, but is 
decoupled from those that have heated the dust observable at 
8 jum or from whose emission is detected in the FUV. 

With a NIR image obtained with the Canada-France- 
Hawaii Telescope, we have identified the presence of a weak 
bar in NGC 3147, a galaxy classified as non-barred in the opti- 
cal. This stellar bar acting on the gas produces gravity torques. 
A portion of the gas present in the inner '^CO(l-O) ring is trail- 
ing the bar, while part of the gas leads it. The gravity torques 
are negative inside a radius of ~2 kpc, while they become zero 
or positive outside. In the outer '^CO(l-O) spiral the torques 
almost cancel out on average, but are positive. In the more 



clumpy inner '^CO(2-l) ring, the predominant torques are neg- 
ative. We interpret these results by identifying the resonances 
with the bar. The inner CO ring would correspond to the ultra- 
harmonic resonance (r = 2 kpc) just inside corotation, and the 
outer CO spiral to the spiral structure emerging from the bar, 
beyond its corotation, which we estimate is located at ~3 kpc. 
In the inner CO ring, the gas is inside corotation due to the neg- 
ative torques there; it loses angular momentum and is flowing 
inwards. 

NGC 3147 is not the first case in the NUGA sample where 
inflowing gas h as been found: NGC 27 82 (Huntetal. 20081), 
and NGC 6574 dLindt-Krieg et al.l2008l) also show this feature. 
However, only NGC 2782 has a strong and significant inflow, 
due to a nuclear bar embedded in the primary one and to gas 
aligned with the nuclear bar The amount of inflowing gas in 
NGC 3147, quantified from the gravity torques computation, is 
modest but could still be sufficient to feed the Seyfert 2 nucleus 
at a rate of about 0.5 Mgyr"'. Although higher spatial resolu- 
tion is needed to determine whether the gas is actually feeding 
the AGN, we are observing the gas inflow to the ~ 100-200 pc 
central region. 
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